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Abstract. The influence of tensor forces on hyperfine splitting and decay widths of 
mesons are investigated within the limits of quasirelativistic potential quark model. As 



the starting the system of the Rarita-Schwinger equations is used. We use a potential 
motivated by QCD with a mixed Lorentz structure. 

1. Introduction 

As well known, the heavy meson spectrum can be successfully obtained starting from 
a potential model [1-3]. The describing of the fine and hyperfine structure including decay 
widths is problematic in framework of unified approach. Many authors with different types 
potentials for quark- ant iquark interaction describe only separate properties of mesons. 
Following papers [4-6] we analyse the Lorentz structure of the potential with reference to 
uu system, charmonium and bottonium hyperfine splitting. It is well known that tensor 
forces bring about mix wave with different angular momentum [7]. Majority of authors 
consider that contribution of tensor forces in spectrum of bound states of quarkonia is 
negligible [8,9]. But, it was shown [3] that D waves give 1 - 5 % contribution for energy 



spectrum in ground triplet state. 

In the present paper we will use model which is based on the quasirelativistic Breit- 
Fermi Hamiltonian [10] and we will be analyze an influence of tensor forces with reference 
to two-quark system. Singlet states are described by Schrodinger equation and for triplet 
states we use the system of Rarita-Schwinger equations [7]. 

2. Hyperfine splitting 

It is widely accepted that the interaction between two quarks (or heavy quark - 
antiquark) consist of a short range part describing the one-gluon-exchange and a infinitely 
rising long-range part responsible for confinement of the quarks. This also gives good 
approximation to the lattice potential. We use 

V = V v {r) + V s (r), (1) 

where 

Vv = ~^-, V s (r) = kr. (2) 

(The color factor 4/3 has been absorbed in the definition of a). Wilson loop techniques 
suggest that the confining potential should be taken purely scalar, but relativistic potential 
calculations which have been published [11-13] show a need for (some) vector confinement. 
Maintaining Vv(r)+Vs(r) unchanged we do allow a fraction of vector confinement [4,6,13]. 
Thus 

VV(0 = -— + A/r, V s {r) = f3 s r {f3 v + f3 s = k). (3) 



Then 



V = (-^ + for) + fa. (4) 
r 



The confining potential transforms as the Lorentz scalar and vector potential transforms 
as the time component of a four-vector potential. As we can see, the choice of Lorentz 
structure of potential for quark-antquark interaction is important model for study of spin 
effects [13-17]. 

We consider the Breit-Fermi Hamiltonian (in case of equal masses m\ = m 2 = m ): 

H = — + V + H LS + H ss + H T (5) 
m 



with spin-orbit term 



S = Si + S2 is the total spin of bound state and L = r x p is the relative orbital angular 
momentum of its constituents, the spin-spin term 

H ss = -^AV v (r)S 1 .S 2 , (7) 



and the tensor term 



^ 1 AdVy d 2 Vy 

12m 2 r dr dr 2 ' 



where 

5 12 = 12[(5!.n)(5 2 .n) - (5i.S 2 )/3] = 3(<n.n)(cr 2 .n) - (<n.<7 2 ), (9) 

n = r/r, 5* = tr f /2. 

For bound-state constituents of spin Si = 5 2 — 1/2, the scalar product of their spin, 
S1.S2, is given by 



S1.S2 



— I for spin singlets, 5 = 0, 
+ \ for spin triplets, 5=1. 



Taking picture (f|), then (|6|), (|7|) and (|j) yield 

^ = ^ + do) 

2m z r 

tf 5S = -^[47ra5(r) + ^(fl^), (11) 

6m z r 

T 12m 2 ^ r 3 r ^ 12 ^ 
Two-quark system may exist in the spin singlet and the triplet states (Table 1). 
If 5 = 0, then necessarily L = J and eigenfunction is of the form (/(r)/r)3^joj, 
where 

iVjls = E( L WI^M)(y im (fi)|5/i). (13) 

mfi 

Since 

5|00 >= 0, 

5 12 = 2[3(5.n) 2 -5 2 ] (14) 
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we get 

S 12 yf LO = S 12 Y Jm (Q)\00>= 0. 
Consequently, f(r) satisfies radial equation 

g^-^-^j/.a (15) 

where l U c = m 1 ^; 1 V c = V v (r) + V s (r) - (3/4)V S s and k 2 = mE. If 5 = 1 and P = 
(— 1) J+1 , the only possible values of L are J + 1 (unless J = in which case there is just 
the one value L = 1). 

The wave function for ground triplet state of qq system with negative parity (P = — 1) 
is a mixture of state 3 Si and Z D\ and may be put in the form 

4 = iP s + iJ> D = - r u(r)y{ m + ^w(r)yl 21 . (16) 

The input D-waves to the wave function of vector meson makes part of percent. But the 
input D-wave to the spectrum of energy makes some percent. It is explained by including 
of the interference member to the energy spectra. Knowing the wave function it is possible 
to calculate input of every component to the energy level. We have 

E = (^ s + tfj D \H\^s + ^ D ). 

Then the energy value can be split into components 

E = E s + E D + E$Di 

where Eg, Ed and represent inputs of S, D-wave and of the interference wave, respectively. 
Now, for S = 1 state we have the wave function fll6|) we know that, 



5l2Xl0 = v8^21 

and 

5i 2 J ; i 1 2i = Vo^oi ~ ^121- 

Futher 

x 2 ^ 01 = o, L 2 y l l21 = eyl 21 . 
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Since 



H 



1 d 2 



+ ^+ 3 V c + V LS L.S + V T S 12 , 



m dr 2 mr 2 

and S1.S2 = 1/4; L.S = —3, S12 = —2 for S — 1. Then the equation (ff — E)tp = is 
equivalent to the Rarita-Schwinger system 

, 1 d 2 



m c?r 2 



E + 3 1/ c 1m + v^Vrw = 0, 



[--Ti - ^ + + 3 K - 2^ T - 3^ L5 ] W + \/8V T u = 0, 



m <ir 2 



where 



3 K - V v + V s + ~V S s- 
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The system ([IT]) we rewrite in the matrix form 



it 



with 



/i = h + W 



(17) 



and 



/ 1 d 



— r^r h/cr-F 

m ar r 

1 \/2/3v 



1 A 



1 V2/3 



v 



m 2 r 3 



m 2 r 6 



m 2 r 6 

1 d 2 6a, 1 
H 7; Vkr 
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m dr 2 mr 2 r 



W 



§AV T 



where 



§AV T \AV S S - j^AVt - 3AV LS ) 

3a . _ 3a 



m z r 3 2 



ma 



ss 



5(f), AVts 



3m 2 ' 2m 2 r 3 ' 
Then matrix elements of correction terms are 



AV T 



m 2 r 3 



AEr, 



1 



,2na 



+ 



u, 



U-r. 



' 3m s 



5{r))u n dr + w 



a 



V2 



2m 2 r 3 



)u n dr 



a 



V2 



'2m 2 r 3 



2ixa 



'3m 2 



5a 

m 2 r 3 



)w n dr. 



(18) 
(19) 



(20) 



We put a = \a s . The QCD running constant a s is determined according to the formula 
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a s (q 2 ) = 12tt/[(33 - 2 nj )ln(q 2 /A 2 )}, 

where A = Aqcd = l^OMeV; rij = 3 for light and mixed mesons; rij = 4 for cc- and bb- 
quarkonium. For Cornell potential hyperfine splitting of mesons may be described with 
a rather accuracy if q — 2/i (/i - reduced mass) is chosen [18]. 

In fact AE mn have singularities (5(r)) and r~ 3 . Therefore one has to introduce smeared 
S function and cut of for r~ 3 to weaken singularities. Our calculation were made with 
tensor forces and without tensor forces. We have used the following parameters of poten- 
tials: (3 V = 0.001GeV 2 ,(3 s = 0.179GeV 2 for mz-systems; (3 V = O.OAGeV 2 , (3 S = O.UGeV 2 
for charmonium and bottonium, a s (cc) = 0.38, a s (bb) = 0.24, a s (uu) = 0.54. The quark 
masses is: m c = lAGeV,m b = A.7GeV and m u = 0.33GeV. In the works [13,16] it was 
pointed out that the best agreement is obtained when (3 = 0.18GeV 2 . During the de- 
scription of gg-system we were changing the parameter (3y (when /3y + Ps = 0.18GeV 2 ) 
to achieve the agreement with experimental mass splitting of 1 -states. Tables 2-4 list 
calculated results. We also display experimental magnitudes [19]. 

Tables 2-4 present the mass spectra of pseudoscalar and vector mesons. They also 
give the size of D-wave admixture in the wave function of vector meson and the input of 
interferential member E$d- 

We also have calculated root-mean-square radii and results are provided in Table 5. 

Squared mean value radius 0.7 fm meet the requirements of quark- ant iquark pair cre- 
ation (uu) (string break). It is necessary modify the potential model taking into account 
the opening of a new channel for those conditions where the given value is overlapped. 

3. Leptonic decay of heavy quarkonia 

For the leptonic decay widths of two-quark system we shall consider decay of 3 S (vec- 
tor) states into e + e~ pairs. The leptonic decay width of system Mqq — * e + e~ is calculated 
from the Van Royen Weisskopf formula [22] 

f( 3 ^i-e^) = %^|i?(0)| 2 . (21) 

Where M qq - is mass of vector meson, Q is quark charge, a em is the fine structure constant 
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and -R(O) is the radial wave function at the origin. 

The formula ([H]) is based on the notion that constituent quark-antiquark pair anni- 
hilates into a single virtual foton, which subsequently gives rise to a leptonic pair. 

The decay widths of the vector bb and cc quarkonia into charged lepton [23] 

r(^ c + c -) = f(i-152iK>). (22) 

As Eichten and Quigg have pointed out [24] the QCD correction reduces the magnitude 
of T significantly, however the amount of reduction is somewhat uncertain . For vector 
mesons containing light quarks this formula leads to paradoxes [25] . 

In paper [2] Motyka and Zalewski have got formula 

Ty^+e- = F{q)f^\Rm\ (23) 

with F(c) = 4.73.10 -5 for charmonium and F(b) = 2.33. 1CT 5 for bottonium. We have 
calculated decay widths by Van Royen-Weisskopf ( pTD and by (p3|). Table 6 lists results of 
calculation. The calculations of widths were done with tensor forces and without tensor 



forces. Value of the widths which were calculated by formula (|23T ) are given in parenthe- 
ses. 

4. Conclusion 

The analysis of our results presented in Table 2-4 shows that percentage differences 
between the theoretical calculations of mass spectrum for heavy quarkonia and the ex- 
perimental results are 1 - 4 %. Therefore we expect relativistic correction in the range 1 
- 16 % for charmonium, and up to 1 % for bottonium spectrum. Also we have been able 
to describe hyperfine splitting of cc- and 66-quarkonium. For describing mass spectrum 
of light mesons it is necessary to use relativistic potential model. We have calculated 
hyperfine splitting uu system, too. 

The leptonic decay widths suggest that for charmonium theoretical widths, which were 
calculated by Van Royen-Weisskopf formula are systematically higher than experimental 
data. But for bottonium we have obtained lower values. For J/^-meson better decay 
widths are obtained by means of formula (23), which takes into account QCD correction. 



Here the influence of D-w&ves ranges from 4 % for ground state up to 25 % for the second 
excited state, but for values calculated by Van Royen-Weisskopf formula it is from 8 % up 
to 50 %. For T-meson it is opposite. More exact is calculation done by (|2l|), but QCD 
correction reduced the results more apparent. Besides, D-waves contribute less than for 
charmonium: from 4 % for ground state up to 11 % for second excited state. 

The results obtained show that contribution of the D-waves is impossible to neglect 
for considered leptonic decay width of quarkonia. 
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Table 1: Conditions of two fermion system 





Singlet condition 


Triplet condition 




(S = 0) 


(S- 


= 1) 


J/p 


+ 




+ 









'So 






1 








3 Si + 3 D 1 


2 




l D 2 


3 P 2 + 3 F 2 


3 D 2 



Table 2: Hyperfine splitting for the charmonium 



State 


S-wave 


SD-w&ves 


[20] 


[19] 


Esd, 


Pd, 




E theor , MeV 


E theor , MeV 


E theor: MeV 


E exp , MeV 


% 


% 


l l S 


2980 






2980 








3153 


3097 




3097 


16 


0.05 


l 3 Sl - 1% 


173 


117 


110 


117 






2% 


3642 






3590 






2 3 ,S 1 


3759 


3734 




3685 


3 


0.8 


2 3 Si - 2% 


117 


92 


67 


95 






3% 


4107 












3^ 


4208 


4192 




4040 


1 


1.3 


3 3 5*i — 3^5*0 


101 


85 
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Table 3: Hyperfine splitting for the bottonium 



State 


S-wave 


SD-w&ves 


[20] 


[19] 


Esd, 


Pd, 




Etheor, MeV 


Etheor, MeV 


Efheor: MeV 


E exp ,MeV 


% 


% 




9415 






- 








9462 


9460 




9460 


2.1 


0.004 


l 3 S l -l l S 


47 


45 


46 








2% 


9883 












2 3 S 1 


9911 


9911 




10023 


0.2 


0.04 


2 3 S l - 2 1 S 


28 


28 


26 








&S 


10201 














10224 


10224 




10355 


0.1 


0.1 


3 3 Si - 3% 


23 


23 











Table 4: Hyperfine splitting for (ww) -systems 



State 


S-wave 


SD-waves 


[20] 


[19] 


Esd, 


Pd, 




E theor , MeV 


E theor , MeV 


E theor , MeV 


E exp , MeV 


% 


% 


I 1 S 


140 






140 








674 


640 




770 


4 


0.02 


1 3 S 1 - 1% 


534 


500 


923 


630 






2% 


1134 






1300 






2 3 5 1 


1564 


1543 




1450 


1 


0.04 


2^! - 2% 


430 


409 


411 


150 
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Table 5: Squared mean value radius 



nL 


[21] 
(r 2 ) 1 / 2 , fm 


Our results 

(r 2 c y/ 2 ,im 


[21] 

(rf) 1 / 2 , fm 


Our results 
(r 2 ) 1 / 2 , fm 


1 S 


0.43 


0.433 


0.24 


0.256 


2 S 


0.85 


0.847 


0.51 


0.552 


3 S 


1.18 


1.182 


0.73 


0.768 



Table 6: The leptonic decay widths of heavy mesons 



State 


SD-waves 


S-wave 


[2] 


[26] 


[27] 


[19] 




T theor ,,keV 


T theor ,,keV 


T theor ,,keV 


T theor ,,keV 


T theor ,,keV 


Texp. ,ke V 


J/iplS 


7.8 (5.41) 


8.2 (5.63) 


4.5 


4.24 


8.0 


5.26 ±0.37 


ifj'2S 


3.7 (2.59) 


4.0(2.79) 


1.9 


1.81 


3.7 


2.12 ±0.18 


*fj"3S 


2.6 (1.82) 


2.9 (2.01) 




1.22 




0.75 ±0.15 


T1S 


1.14 (0.96) 


1.20 (1.01) 


1.36 


0.85 


1.7 


1.32 ±0.04 


V2S 


0.58 (0.49) 


0.63 (0.53) 


0.59 


0.38 


0.8 


0.52 ±0.03 


T"3S 


0.44 (0.37) 


0.49 (0.42) 


0.40 


0.27 


0.6 


0.48 ±0.08 
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